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Abstract 



We formulate the double exchange (DE) model for systems with tig orbital 
degeneracy, relevant for hole-doped cubic vanadates. In the relevant regime 
of strong on-site Coulomb repulsion U we solve the model using two distinct 
mean-field approximations: Hartree-Fock, and mean-field applied to the for- 
mulation of the model in the Kotliar-Ruckenstein slave boson represenation. 
fj ' We show how, due to the relative weakness of the DE mechanism via tig 

degenerate orbitals, the anisotropic C-type antiferromagnetic metallic phase 
and the orbital liquid state can be stabilized. This is contrasted with the DE 
mechanism via e g degenerate orbitals which stabilizes the ferromagnetic order 
in the hole-doped regime of manganites. The obtained results are in strik- 
ing agreement with the observed magnetic structures and the collapse of the 
Z^. orbital order in the doped Lai-^Sr^VOs, PTi- x Ca, x Y03 and Ndi-^Sr^VOa. 

o 

\Q ■ 1 Introduction 

p' 

Recently orbital degrees of freedom in the strongly correlated electron systems have 
attracted much attention pQ due to their crucial role in the stability of the various 
magnetic phases in the Mott insulators P] or in the partial explanation of the 
Mott metal-insulator transition itself pj]. Typical examples of such systems are 
the transition-metal oxides with partly filled e g or t^g degenerate orbitals. The 
remarkable feature of these oxides is the large on-site Coulomb repulsion U 3> t, with 
O ■ t being the effective electron hopping. In case of the stoichiometric (undoped) oxides 

the Coulomb repulsion between electrons supresses charge fluctuations and leads to 
the effective low-energy ex J=At /U superexchange (SE) interactions between spin 
and orbital degrees of freedom. In addition, the atomic Hund's interaction Jh <£/ 
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C$ • aligns the spins of the electrons occupying the degenerate and almost degenerate 

orbitals on the same site and should be taken into account in the realistic SE-type 
models. Interestingly, it is the Hund's interaction which to large extent stabilizes 
particular spin (magnetic) and orbital order. It should be stressed that though the 
e g orbitals usually order as in e.g. LaMnOs |2J, it is rather not the case of the 
t2 g orbitals which can form a disordered quantum orbital liquid (OL) as in e.g. 



LaVC>3 Hj. The different kinds of the orbital state are concomitant with various 
anisotropic magnetic phases such as ^4-AF [ferromagnetic (FM) planes staggered 
antiferromagnetically (AF) in the other direction] in LaMn03 [2] or C-AF [FM 
chains staggered AF in the other two directions] in LaVC>3 |3]. 

Doping with holes destroys the insulating state, modifies the SE interaction 
and hence can also modify the magnetic order and orbital state stabilized in the 
undoped case. In such systems the motion of holes is strongly affected by: (i) 
intersite SE interaction ex J, (ii) on-site Hund's interaction oc Jh, which is captured 
by the Kondo-lattice model. In the realistic limit in the transition metal oxides 
Jh ^* t 3> J this model reduces to the double exchange model (DE) Qp. However, 
for the above discussed class of the transition- metal oxides, orbital degeneracy needs 
also to be taken into account in such a model. This was done in the seminal work 
of van den Brink and Khomskii 6] for electrons with the e g degeneracy, leading 
to a naively counterintuitive picture of the DE mechanism. One can thus pose a 
question how the DE mechanism would be modified in the case of the t% g orbital 
degeneracy, as the SE interactions for the undoped oxides differ qualitatively for the 
tig and e g cases. Hence, in this paper we want to answer three questions: (i) what 
is the nature of the magnetic order stabilized by the DE mechanism for correlated 
electrons with t 2g orbital degeneracy, (ii) what is the nature of the orbital order, 
and (iii) how do these results differ from the e g and the nondegenerate case. 

The paper is organized as follows. In the following chapter we introduce the 
DE Hamiltonian with ti g orbital degrees of freedom. Then we solve the model 
using two different mean-field approximations and we show that the metallic state 
coexists with the C-AF order and OL state for the broad range of hole-doping. 
Next we discuss the results, i.e.: the validity of the approximations, the generic role 
of the tig orbitals in the stability of the above phases, the physical relevance of the 
model by comparison with the experiment. The paper is concluded by stressing: 
the distinct features of the DE mechanism via e g and via t 2g degenerate orbitals, 
and the crucial role of the Coulomb repulsion U. 

2 Realistic double exchange Hamiltonian 

We start with the realistic semiclassical DE Hamiltonian with t 2g orbital degrees of 
freedom, relevant for the hole-doped cubic vanadates 0: 

hde = -u y y A«j - 2 ^ - u * Y && - Y *$j +J Y. s ^ s r co 

where: the hopping amplitude t=l; the restricted fermion creation operators a\ — 
a I(l — "v>i)i b\ = b\(\ — a\a,i), where aJ(&T) creates a spinless electron at site i in 
yz(zx) orbital; and Si are core spin S = 1/2 operators of t 2g electrons in occupied 
xy orbitals. We introduce the u x and u y variational parameters to be defined later. 
The Hamiltonian has the following features: (i) the first four terms describe the 
kinetic energy of the electrons in the degenerate yz and zx orbitals which can hop 
only in the allowed (y, z) or (x, z) plane to the nearest neighbour (nn) site i providing 



there are no other electrons at site i in these orbitals (U — oo assumed implicitly), 
(ii) the last term describes the AF coupling between core spins at nn sites in the 
(x, y) plane due to the SE interaction originating from the electrons in the always 
occupied xy orbitals (n xy = 1 during doping), (iii) the SE interactions due to the 
itinerant electrons in yz and zx orbitals are neglected. 

The variational parameters are defined as u x — coa(9 x /2) and u y — cos(9 y /2), 
where 9 X (9 V ) is the relative angle between core spins in the x (y) direction. This 
follows from the Hund's rule which aligns the spins of the itinerant electrons with 
the core spins and which is not explicitly written in Eq. QJ but enters via u x and 
u y parameters. Then the last term of Eq. (QJ can be written as: 

Ese = J J2 S J -S, = ^ 3 (^ + u 2_ 1)i (2) 



(y'JPifi 



where L 3 is the number of sites in the crystal. Eq. Q together with Eq. @ shows the 
competition between the DE mechanism allowing for the hopping in the underlying 
FM background with the SE interaction supporting the AF order. 

3 Numerical results 

The ground state of Eq. Q was found using two distinct mean-field approximations: 
Hartree-Fock, and mean-field applied to the formulation of the model Eq. (QJ in 
the Kotliar-Ruckenstein slave boson representation [H]. 

3.1 Hartree-Fock approximation (HFA) 

To facilitate the treatment of the hole doped electronic systems we enlarge the 
original Fock space, which is relevant for the Hamiltonian Eq. (QJ, by adding a 
boson creation (annhilation) operator e|(ej), which creates (annihilates) a charged 
hole. Then, in order to get rid of the nonphysical states in the enlarged Fock space 
we need a local constraint: 

Vi a\ai + b\bi + e\ei = l. (3) 

Note that the boson controlling double occupancies is not introduced since such 
configurations are forbidden for U — oo. In order to solve Eq. (QJ with the constraint 
Eq. (jHJ in a simple way, we introduce mean-field approximation for the above 
constraint and replace the following combination of fermion operators by averages: 

Vj b\bi -> (b\bi) =: 5 and < <5 < 1, (4) 

Vi a\ai — > (a-lai) = 1 — x — 8 and < 5 + X < 1, (5) 

where in Eq. 10 we use constraint Eq. Q), where (eTe,) = x and x is the number of 
doped holes per site. Besides, the inequalities in Eq. <@} and JU follow from the fact 
that we are discussing the hole-doping regime, i.e. with not more than one itinerant 
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Figure 1: Results of the minimisation of Eq. Q (with the constraint Eq. 10): (a) 
Magnetic phase diagram at T — 0. Under increasing hole doping x, the C-AF order 
changes via the canted Ca** phase, with spins canted on two bonds, to a homogenous FM 
phase. The grey horizontal line depicts the realistic value of SE J = 0.18 in LaVOs; (b) 
Competing: correlated kinetic energy (-Ekinetic), and magnetic energy (Use) as a function 
of hole doping x for the fixed value of SE J = 0.18. Kinetic energy for uncorrelated 
electrons (-Eband) and the total energy of the system (.Etotal) is also shown. 

electron per site. We call this mean-field approximation an HFA approximation 
since it is equivalent to decoupling the interaction term in an HFA way: 

(U-3J H ) 4aib\bi ~ (17-3 J H ) ({a\<H)b\bi+4^\bi)-{aia\){b\h^ 



(6) 

and then making the U — ► oo limit [where U — 3Jh is the energy of the t\ g high spin 
state]. Note that the above used HFA does not allow for the alternating orbital (AO) 
order but we do not expect such an order to persist during doping. This is because 
the kinetic energy of the electrons in the AO is rigorously zero for U — oo and the 
vanishing hopping between different orbitals (Hj. Then Hbe can be diagonalized 
and the total energy of the system in the mean-field approximation is: 



E H fa/L=(1-S) 



£(ky,k z ;uy) + (x+8)' 
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£{k x ,k z ;u x ) + E; 



SE 



(7) 



where the kinetic energy renormalizing factors follow from the HFA. The sums go 
over the mesh of k's belonging to the Fermi volume Sf which "accommodates" 
1 — x — 8 or 6 electrons per site (respectively) and whose shape depends also on the 
band dispersion relations for the tight-binding model renormalized by the underlying 
magnetic order: 

£(k y , k z ]u y ) = —2{u y cosk y +cosk z ) and £{k Xl k z ;u x ) = —2(u x cosk x +cosk z ). (8) 

Since we are calculating separately the energies of the yz and zx bands we need an 
additional constraint for the Fermi levels £f as the two bands fill up gradually: 



£f(1 - x - 5, Uy) = £f(8, u x ) 



(9) 



For fixed values of SE J and hole doping x we minimise numerically the to- 
tal energy Eq. (with the constraint Eq. JUJ) with respect to the u x and u y 
variational parameters, which determine the magnetic structure. The resulting 
magnetic phase diagram at T = is shown in Fig. ^a). In Fig. QIb) we show the 
competing kinetic and magnetic (SE) energies as a function of hole doping x for 
the realistic value of SE J = 0.18 and obtained in the minimum of the total energy. 
Let us also note that for all of the discussed values of J and x it was found that 
gmin _ ( c ). Ci j — (d\di) = (1 — x)/2, i.e. the number of electrons in both orbitals is 
the same and the ferrorbital order (FO) is not stabilized. 

3.2 Kotliar-Rukenstein slave boson approach (KRA) 

We rewrite the Hamiltonian Eq. (QJ using Kotliar-Ruckenstein slave boson represen- 
tation [HI adopted to the orbital case. So we enlarge the Fock space by introducing 
three auxiliary boson fields and decouple the electron creation operator a\(b\) into 
a fermion creation operator f}(g{) carrying the orbital degree of freedom (orbital 
flavour), a boson creation operator p](ql), and a boson annihilation operator e, 
carrying the charge degree of freedom (i.e. ej creates a charged hole): 

4 = fUt h\ = g\yi (10) 

where 

t Pi^i J t Qi^-i / 1 1 \ 

»i'= I =. and vl-= , ■ =■ (ii) 

y (1 - e^ei-qlqi)(l-piPi) y (1 - e'<e l -plp l ){l-q]q l ) 

Applying this to the Hamiltonian Hue we get: 

Hde = - % J] fjfjzjzj - Yl flfi4 z J + ^SE 

i,j\\v <.i||i 

- u x Yl dldivhi - Y dtmtyj, (12) 

where we also need constraints to get rid of the nonphysical states in the enlarged 
Fock space: 

V; PiPi+qjqi + eje; = 1, V, ptp i = ftf i and V 4 qjqi=gtgi- (13) 

Next we introduce mean-field approximation: 

(pl^elpi) _ x 



Vj z\zi -> (z\zi) = 



(1 - eJe t - q\qi){l - p\pi) %a 



w t / t \ k f e t ejft) a; 

(l-eJe i -plp i )(l-qlq l ) x b 



0.3 



0.2 



0.1 



0.0 



- 


C-AF 


Ca** 


(a } 


- J=0.18 






- 


C-AF 








•r . . . . i . . 


... i ... . 


FM 

■ ' ' ■ 





-1.2 



-1.6 



W 

***** **** + : 



00 h « I ^ ^ 

▼ T * ## ^E T<inetic 

-04 ■ '''j'l, Sand T Sotal 

I -0.8: "• : »i:::r 



iii» ij 



ti^' 



0.0 0.1 0.2 0.3 0.4 0.5 
x 



0.0 0.1 0.2 0.3 0.4 0.5 



Figure 2: Results of the minimisation of Eq. 1151 (with the constraint Eq. 10): (a) 
Magnetic phase diagram at T = 0. Under increasing hole doping x, the C-AF order 
changes via the canted Ca** phase, with spins canted on two bonds, to a homogenous 
FM phase. The grey horizontal line depicts the realistic value of SE J = 0.18 in LaV03. 
Note that for x = the obtained magnetic order consists of decoupled AF (x, y) planes 
[after Ref. 0]; (b) Competing: correlated kinetic energy (-E kinetic ), and magnetic energy 
(.Ese) as a function of hole doping x for the fixed value of SE J = 0.18. Kinetic energy for 
uncorrelated electrons (-Eband) and the total energy of the system (.©total) ' s a ' so shown. 

where we use Eq. II 1811 and introduce: (1 — pjpi) — x a the mean number of holes 
per site in orbital yz, (1 — qjqi) — xt, the mean number of holes per site in orbital 
zx, and (ejej) = x the mean number of holes per site. Furthermore, introducing 
parameter S defined in the same manner as in Eq. (|4I5I) yields x a = x + 5 and 
Xb = 1 — S. Then the total energy of the system in the mean-field approximation is: 
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where the electronic bands are defined in the same manner as in Eq. (jHJ and the 
constraint Eq. (jHJ also holds. 

As in the HFA we deduce the resulting magnetic structures by minimising the 
total energy Eq. i[l5ll with the constraint Eq. Jsj. The resulting magnetic phases 
are shown in Fig. |2{a) and the competing energies at fixed SE J as a function of hole 
doping x are shown in Fig. Efb) . We should also stress that the number of electrons 
in both orbitals does not change during doping and is equal to 5 mm = (1 — x)/2. 
Hence, the kinetic energy renormalizing factors are equal to 2x/(x + 1), the well- 
known Gutzwiller factors. 

4 Discussion of the results 

4.1 Why Hartree-Fock approximation fails 

Let us first discuss the validity of the two mean-field approximations. For the HFA 
we infer that for small doping (x — ► 0) the correlated kinetic energy (-Ekinetic in 



Fig. IHb)) does not tend to zero and instead equals 1/4 of the uncorrelated kinetic 
energy for x = 0. This is a remnant of the fact that in the HFA the probability 
for existence of the hole at each site is equal to 25%. Thus, we conclude that the 
configurations with double occupancies are not suppressed and that HFA fails, at 
least in the low hole doping regime. The reason why we cannot use HFA, unlike in 
e.g. the undoped antiferromagnet, is that the orbital state, which is favoured here 
by the DE mechanism, is not an ordered state. If we were to get the FO or AO 
state, HFA would give physically relevant results. 

4.2 Stability of the C-AF: generic role of the degenerate t 2g 
orbitals 

We concentrate on the results obtained within KRA which give the (qualitatively) 
right value of the correlated kinetic energy of the Hamiltonian Eq. (T1J, cf. Fig. Etb). 
The C-AF order and a concomitant one-dimensional (ID) metallic state is stable for 
a broad range of the phase diagram parameters, in particular for the realistic value 
of SE J — 0.18 and < x < 0.28 hole doping. This means that the DE mechanism 
does not win in the (x,y) planes, where the SE mechanism takes over. The first 
explanation of this phenomenon would be to ascribe the role of the driving force for 
the stability of the C-AF order to the Coulomb repulsion U: the correlated kinetic 
energy (^kinetic in Fig. E3[b)) is much smaller than the bare band energy (-Eband 
in Fig. Bib)) and hence the possible kinetic energy gain in the homogeneous FM 
state is much smaller in the correlated regime than in the conventional free electron 
regime. However, for hole-doped manganites U, Jh and t are of the same order as 
for hole-doped vanadates [Uj]. But in the manganites for x > 0.18 hole doping the 
homogeneous FM order, concomitant with the metallic state, is stabilized mainly 
due to the DE mechanism Ql]! The only difference is that in the manganites the 
itinerant electrons hop between degenerate e g orbitals and in the vanadates between 
degenerate t 2g orbitals. Besides, the AF interaction between core spins is truely 
three-dimensional in manganites and not two-dimensional (2D) as it is the case of 
our model. Thus the relative weakness of the DE mechanism in the Hamiltonian 
Eq. {U can be attributed to the specific features of the t 2g electrons: the strictly 
2D and strictly flavor conserving hopping between these degenerate orbitals 0. 
This generic role of the t 2g orbitals in the stability of the C-AF order is threefold: 
(i) The AF SE interactions are rigorously 2D because for xy electrons the virtual 
hopping processes leading to the SE interactions can happen only in the (x,y) 
plane. Hence it makes possible: the hopping in the z direction without destruction 
of the magnetic order, and thus the gain of the kinetic energy without any loss in 
the magnetic energy. Though this is rather a specific feature of our simplified SE 
interactions than of the DE mechanism itself it should be stressed that the physical 
setting in which we can have the DE via t 2g degenerate orbitals is such that the AF 
core spin interactions are 2D. 

(ii) The degeneracy of the orbitals reduces the correlated kinetic energy and 
hence makes it energetically unfavourable to get homogenous FM order. Cf. Fig. 
2(a) of Ref. (3 where for the system without orbital degeneracy with holes doped 



only into zx orbital, i.e. with the Hamiltonian: 

U zx = -u x 53 4^ - 5] b\b 3 + ^L 3 (ul +u 2 y -l), (16) 

njp *>iP 

we obtain that for the realistic value of J = 0.18 the FM order is already stable 
for hole doping x > 0.1. Note that this comparison is valid only if the degenerate 
orbitals are in the OL state which indeed is the case, cf. discussion below. 

(hi) Since each of the ti g orbitals has a different inactive axis [i.e. the cubic 
direction in which the electrons cannot hop] this means that if one of these orbitals 
is "engaged" in the AF SE interactions [cf. (i)] than the other two orbitals have 
one active axis in the direction perpendicular to the plane with AF bonds. Hence, 
hopping in the z direction, assumed in point (i), is indeed possible. 

4.3 Stability of the orbital liquid 

The orbital state which is concomitant with the obtained magnetic phases is the 
OL since: (i) we obtained x a = Xb for all of the magnetic phases, (ii) x a = x^ => 
(T?) := (p]pi) — {QiQi) — where T* is the z component of the pseudospin orbital 
operator j^j, (in) {a]bi) = (frja,) = =>- (Tf) — (Tf) —Q for the x and y components 
of the pseudospin orbital operator [§J . The reason why none of the ordered state are 
realized in nature is twofold. Firstly, as already described, the AO state prohibits 
hopping due to the vanishing hopping integrals between different orbitals. This 
mimics the confinement of holes in the Neel AF, but stays in contrast with the 
behaviour of the e g orbitals. Thus we cannot gain any kinetic energy upon doping 
and this state is unphysical. Secondly, the FO state is also unphysical. A priori, 
due to the Pauli Principle the electrons in the polarized state can better avoid 
each other than in the unpolarized one. Hence the correlated kinetic energy of 
the FO state is enhanced in comparison with the OL, cf. Fig. 8 (with 7 = 0) of 
Ref. [Hj. Though, the condition of the equilibrium for the whole system yields that 
the chemical potentials (Fermi energies) of the two fermionic subsystems (one with 
electrons in yz orbitals and the other one in zx orbitals) should be equal. For the 
FO state to be realized in nature it means that the DE mechanism should raise the 
bottom of one band so that the electrons can pour into the other one. However, the 
DE is too weak in our system and this does not happen. 

4.4 Physical relevance of the model: comparison with the 
experiment 

As already pointed out our model describes the physics present in the hole-doped 
cubic vanadates such as Lai-^Sr^VOs where the following phases were observed 
experimentally upon changing the hole doping x in T — > 0: C-AF, AO and insulating 
one for x < 0.178 hole concentration, C-AF and metallic one for 0.178 < x < 0.26, 
and paramagnetic and metallic one for 0.26 < x < 0.327 [12|. Rather similar phases 
were observed experimentally in Pri-^Ca^VOs and Ndi-^Sr^VOs |12j . Hence our 
results stay in well agreement with the experiment, as we did not expect to explain 



by our model Eq. (QJ the filling control metal-insulator transition present in the 
system, but merely the coexistence of the metallic phase, the OL state and the C- 
AF magnetic order in the above mentioned doping range. However, let us note that 
including the SE interactions between the itinerant electrons would stabilize the ID 
OL in the z direction and AO in the (x, y) plane [3]. Thus, if we included the realistic 
SE in the Hamiltonian Eq. (QJ we would get the insulating behaviour for finite hole 
doping due to the AO order in the plane which would act as a string potential 
for holes moving in the z direction. This suggests that further investigation of the 
even more realistic DE model with t 2g orbital degrees of freedom and SE interactions 
between itinerant electrons would not only uncover some new interesting physics but 
also could explain the experimentally observed phases in the doped cubic vanadates. 

5 Summary 

In summary, we have shown the distinct features of the double exchange via t 2g 
degenerate orbitals. These features, attributed to the t 2g orbital degeneracy, are 
the following: 

(i) Similar to the electron-doped e g degenerate orbitals j^] the orbital degeneracy 
naturally allows for the coexistence of the anisotropic magnetic phases and the 
metallic phase. This is a generic feature of the DE mechanism via degenerate 
orbitals. 

(ii) In contrast with the hole-doped e g case ^T] the homogenous FM order is not 
stable for large range of doping in the discussed here hole-doped regime of ti g case. 
This is due to the different geometry of the t 2g orbitals causing the hopping being 
strictly 2D and strictly flavour conserving, hence weakening the DE mechanism. 

(iii) Upon hole-doping the orbitals form the OL state due to the weakness of 
the DE mechanism. Hence, they do not order as it is typical for the electron-doped 
e g systems (j^j. However, it should be stressed that the type of the orbital state for 
hole-doped e g systems is still under discussion [9|. 

Besides, the correlated behaviour of the itinerant electrons plays a crucial role 
in the stability of the obtained magnetic and orbital phases. This is especially 
dramatic in the orbital sector, where the orbitals form the orbital liquid and the 
HFA [suitable either for weakly interacting electrons or for the ordered states] fails 
qualitatively in giving the right values of the kinetic energy. Altogether these are 
striking results which can serve as a basis for further investigation of the correlated 
phenomena in hole-doped systems with t 2g orbital degeneracy. 
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